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Abstract 



In low magnetic field, the stacked, triangular antiferromagnet CsCuCla 



'^ I has a helical structure incommensurate (IC) in the chain direction. The IC 

C I wavenumber (from neutron-diffraction experiments) decreases with increasing 

Q ' field transverse to the chains, as predicted by classical theory, but then it has 

a plateau almost certainly caused by quantum fluctuations. Linear spin- 
wave theory fails because fluctuations have particularly large effects in the IC 
phase. An innovative phenomenological treatment of quantum fluctuations 
O ! yields a plateau at approximately the observed value and the observed fields; 

^^ ' it predicts a transition to the commensurate phase so far not observed. Results 

r^ • depend sensitively on a weak anisotropy. 

I> , 

Q\ : 64.70.Rh, 75.10.Jm, 75.25.+Z, 75.50.Ee 
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I. INTRODUCTION 

Compounds of the ABX3 family {A = Rb, Cs; B = Mn, Fe, Co, Ni, Cu, V; X = CI, Br, I) 
figure prominently in the study of phase transitions in low-dimensional systems. Much of 
the interest in the magnetic-field behavior of compounds like CsCuCls arises because they 
are physical realizations of models related to the triangular antiferromagnet (TAFM). The 
TAFM ground state is both continuously and discretely (two-fold) degenerate, even in the 
presence of a magnetic field H (with magnitude H less than the saturation field Ho,), unlike 
that of the square-lattice AFM. In a field, thermal fluctuationsuQ in classical TAFM models 
and quantum fluctuationsH'H break the continuous degeneracy (which is nontrivial because 
it is not due to a symmetry of the Hamiltonian) in the same way, both selecting for example 
the colinear structure at if ~ iJs/3 • 

The magnetic properties of CsCuCls (with a Neel temperatureQ Tn = 10.7 K) arise from 
the Cu"*"^ ions; to a good approximation, these form a triangular lattice of parallel chains, the 
other ions serving to define the structure. The major interactions, all nearest-neighbor, are 
a ferromagnetic exchange interaction in the chain or c direction, a weaker antiferromagnetic 
exchange interaction between chains (within the a-b planes), and a Dzyaloshinskii-Moriyau 
(DM) interaction also in the c direction. Both exchange interactions are nearly isotropic; the 
latter is frustrated. Recent studiegHM of the structure and of the phase transition giving rise 
to the DM term cite earlier literature on these topics. In the simplified structural model, 
the classical, zero-temperature, zero-field structure is a three-sublattice, ±120° TAFM 
structure in each a-b plane; the spins lie in the nlanes and rotate from plane to plane, 
forming an incommensurate (IC) helical structureEJ. 

Fluctuation effects are likely large in CsCuCls for several reasons: the Cu spin is small 
{S = 1/2), the system is almost one-dimensional (the intrachain interaction is much larger 
than the interchain interaction), the interchain interaction is frustrated, the exchange inter- 
actions are nearly isotropic, and the structure is incommensurate. Experiments in magnetic 
field (difficult because Hs = 30 T) indeed find major effects due to quantum fluctuations. 

CsCuCla in a longitudinal fleld (H || c) appears to be well understood at low temperatures 
T, but the transition at Tn has puzzling featurest3. The discontinuityii^nlJ in the low-T 
magnetization at H ^ 0.4 H^ was shown by Nikuni and ShibaE3E3 to be a novel, fluctuation- 
induced phase transition from the umbrella structure (optimal at small H due to a small, 
easy-plane anisotropytS in the intrachain exchange) to a coolanar structure (optimal at larger 
H due to quantum fluctuations). Further experimentsc3iH^, including neutron-diffractionEll 
and specific-heat measurementsllj near T^ , confirmed their analysis. 

Properties for a transverse field (H ± c) are not well established. Small fields de- 
form the helix, increasing its period; in agreement with experimentc^'cil, classical (mean- 
field) theoryEl'E^ predicts that the IC wavenumber q decreases quadratically as H increases 
and that the curvature increases with T. The structure remains incommensurateEil up to 
H ^ 0.44 i^s! it is unknown at larger fields, where classical theoryE§cJ predicts an inter- 
mediate commensurate (C) phase in which each plane has the same three-sublattice struc- 
ture. Classical theory fails at intermediate fields: near Hs/3 , plateaus are observed in the 
magnetizationEJ m, the ^^^Cs NMR shiftE2l and the wavenumberE^; the ESR measurementsE3 
are not easily interpreted. The plateau in m, as for the TAFMH'El, is due to quantum 
fiuctuationsE3 (the analysis was done for the C state, but the result for the IC state should 



not differ much). Quantum ffuctuations are surely responsible also for the plateau in q, 
but this remains to be demonstrated. Recent specific-heat, magnetization and neutron- 
diffraction measurement sell near Tm suggest maior effects due to thermal ffuctuations: Tm 
increases with field (as in the TAFMlil'Qc3) and a new phase appears. 

The following examines quantum ffuctuations in the incommensurate phase of CsCuCls . 
Section II. A describes the Hamiltonian. Section II. B describes classical results (from solution 
of the Euler-Lagrange equations); surprisingly, the spins do not remain in the a-h planes 
at intermediate fields (as found also by Jensena). Section II. C describes a linear spin- 
wave (LSW) analysis based on the classical results; this fails because LSW theory does 
not find the C-state and IC-state energies to equivalent accuracy. Section III introduces a 
phenomenological treatment of quantum ffuctuations and shows that it works well for the 
C state; the same approach applied to the IC state yields a plateau in the wavenumber 
at approximately the observed value and at approximately the observed fields. Although 
results are rather sensitive to an anisotropy parameter, we conclude that a commensurate 
phase should appear at a field well below the saturation field, likely below 0.5 i^s • 



II. HAMILTONIAN AND ANALYSIS 

A. Hamiltonian 

The Hamiltonian corresponding to the simplified structure is 

'^ = 2^ [ — 2 Jo Sj„ ■ Sj^n+i + 2 ?7 Jo S'j„ S^^n+i — -D z ■ (Sm x Si,„+i) 

in 

+ Jl/Z ^in' ^^n~ 9 I^bH-K- Sin , (1) 

k 

where Sj^ is the spin operator at the z-th site in the n-th a-h plane, z and x are unit vectors 
in the c and a directions, and the k sum is over the six, in-plane, nearest neighbors of the site 
in. The first term (oc Jo) is the isotropic, ferromagnetic exchange interaction between spins 
in nearest-neighbor planes, the second is an anisotropic correction (of easy-plane type) to 
the first, the third (oc D) is the interplane DM interaction, the fourth (ex: Ji) is the isotropic, 
antiferromagnetic exchange interaction between nearest-neighbor spins in the a-h planes, and 
the fifth is the Zeeman energy in a transverse field H. The values of the coefficients have 
been estimated in several articlesQ^MtaJ^; we use Jo = 28 K, r^Jo = 0.24 K, Ji = 4.9 K and 
D = 5 K. We omit the dipole-dipole interaction (which can induce an IC state modulated in 
the planesE3), the anisotropy of the interchain interaction, and several small, related effects 
(namely the displacement of the Cu ions from the c axis, the component of the DM vector 
perpendicular to the c axis, and the z component of the magnetizationgH^cHl) . The saturation 
field iJs = 30 T, above which each spin is aligned with the field, is 18JiS'/((7/iB) in terms of 
the model; the reduced field is h = H/Hg . A longitudinal field (H || c) maintains the axial 
symmetry of zero field, and the DM term can be eliminatedE^O. A transverse field (H ± c) 
breaks the symmetry, and the DM term comes into full play. 

The intrachain exchange term (Jo) favors a state with spins parallel in adjacent layers 
while the smaller DM term favors a rotation by 7r/2 per layer; the result at low fields is 
a helical structure, incommensurate in the c direction. The IC wavenumber q is small 



{^ 27r/71) at h = where the continuous degeneracy of the classical TAFM ground state 
corresponds to a mere shift in the origin of the coordinate system; the classical degeneracies 
remain for a transverse field < H^ , even if the spins are confined to the planes. 

Both the DM term (Inactive for the possible C state in transverse field) and the much 
smaller anisotropy termt3 (responsible for the phase transition in longitudinal fieldt^'tll) favor 
spins in the a-b planes. For the IC state in transverse field, it is natural to assume that the 
DM term confines the spins to the planes at all fields (as it does at small H), and to ignore 
the 7] term. Surprisingly, in-plane spins are unstable at intermediate fields for isotropic 
exchange (as found also by JensenE^); the instabihty persists, but is much reduced, for the 
experimental anisotropy. 

We assume that the three-sublattice structure is maintainedEil'CJ at all fields < Hs . 
A continuum approximationc3 for the c-direction dependence is valid, but we require the 
discrete model because we approximate IC states as high-order commensurate states; in 
principle, we miss some features of IC states, but on the other hand we find that the states 
are not pinned. We use periodic boundary conditions Sj^i+l = ^ji , where j = 1,2,3 is the 
sublattice index and / = 1, ■ ■ ■ , L is the layer index. 



B. Classical analysis 

In the classical approximation, the spin operators Sj„ become classical vectors of length 
1/2, and the Hamiltonian of Eq.(l) becomes the energy function Ec\{{Sin})- With Lagrange 
\ji for the constraints S^; • S^; 



multipliers Xji for the constraints S,-; • Sji — S"^ = 0, the Euler-Lagrange equations are 



- 2 Jo (S,, i_i + S,,+i) + 2r]Jo {S\f_, + S^^%,) z - D (S,- ,+i - S,- ,_i) x z 

+ 6 Ji (S,_i,i + S,+i,0 - 18 Ji Shi -2 \,i Sji = ; (2) 

the solutions are not pinned and so a condition such as y ■ Sn = is necessary. Ref. ^ 
studied the continuum version of Eq.(2) (without the anisotropy term) for in-plane spins. 

The Euler-Lagrange equations have many solutions. The one optimal at a given field is 
too often neither simple nor obvious, and so it was necessary to generate scores of solutions 
and follow them as functions of the field (in steps Ah = 0.01), minimizing the energy of 
each with respect to L ; composite solutions (with increased winding number) gave greater 
accuracy when such was desired. Random starting configurations were generated at various 
values of h and L, relaxed by conjugate-gradient minimization of the energy, and then used 
to start solution of the Euler-Lagrange equations. 

Classical theory cannot explain the plateau in q , but it is a necessary preliminary to the 
LSW analysis and it sheds light on two questions: (1) does a commensurate state intervene 
between the incommensurate and aligned states, and (2) do the spins lie in the a-b planes 
at all fields? The answers depend crucially on the size of the anisotropy parameter 77. 

For isotropic exchange (77 = 0), the optimal structure is incommensurate up to Hs . 
For h < 0.38, there are many solutions (as in the continuum approximationEa), but they 
are well separated in energy, and the spins remain in the planes for the optimal solution 
(the 111 solution of Ref. EH), which evolves continuously from the solution at /i = . For 



h ^ 0.38, the spins break out of the planes (in the optimal solution); this is surprising, 
for each intrachain term in the energy is optimized or neutral if the spins lie in the planes. 



and so is the sum of the interchain and field terms. In contrast, we verified that the spins 
remain in the planes for all fields in the case of a modelEj with a ferromagnetic chain-chain 
interaction. The IC phase for h ^ 0.38 is, however, very complicated: the solutions are far 
more numerous, they differ only slightly in energy for the most part, and the low-energy 
solutions are not the simplest; ground-state energy crossings make it impossible to identify 
with confidence the optimal solution in most of this range. In short, the IC phase is glassy for 
fields in much of the range 0.38 ^ /i < 1, the ground state changing often and unpredictably 
with field, from one complicated configuration to another. 

The weak, easy-plane anisotropytJ r] = 8.6 x 10~^ changes the phase diagram sub- 
stantially, opening a large window for the C phase and destroying the glassy phase. For 
< h < 0.41 , the in-plane IC solution of the previous paragraph is optimal; it extends to 
larger h because of the anisotropy. For 0.42 < h < 0.50 , a simple, out-of-plane, IC solution 
is optimal (but other IC solutions are close in energy); that is, both the DM and anisotropy 
terms acting in concert fail to confine the spins to the planes at all fields. For 0.51 < h < 1 , 
the commensurate (C) state is optimal. All transitions are second-order. If t] were only 
slightly larger (say a few %), then likely the spins would remain in the planes at all fields, 
and the IC phase make a second-order transition to the C phasec3 at h ^ 0.47. 



C. Linear spin wave analysis 

The following describes a linear-spin-wave (LSW) analysis of quantum fiuctuations in 
incommensurate (IC) states. To our knowledge, fiuctuations have previously been studied 
only in IC states which are sinusoidal or nearly so. The analysis was restricted to in-plane 
spins and isotropic exchange (?7 = 0), and so results were obtained only for h < 0.38 . 

A Holstein-Primakoff transformations to boson operators (using a local coordinate 
systemca) and an expansion about the classical solutions give the Hamiltonian as 

n = E,i + ni + n2 + o{^^) ; (3) 

the classical spins are not colinear, and so the expansion parameter is l/v>5 (rather than 
1/S), as in Refs. ^ and ^ The classical angles (pji between the spins and the field are 



determined by minimizing the classical energy E^ (of order S"^) of the A^ spins: 

JV "^ . _^. 
•^^ i=i 1=1 

+ 6 Ji cos{(f)j+i^i-4>ji) -ISJih cos 4>ji . (4) 

The term Tii (of order S^^'^) is linear in the boson operators; it vanishes for the classical 
angles (pji and so is omitted. The LSW Hamiltonian 7^2 (of order S) has the standard form 
(quadratic in the boson operators): 

^2 = -f ^EC^.^ + f EEElA^.'Kk) [6],(k)6^.,,(k) + 6^,(-k)4^,(-k) 



2 3^ ji 2 ^ -i .,^, 



+B,iM^) [h],{\,) 4,(-k) + 6^,(-k) h.yM] } . (5) 
5 



The k sum runs over the N/{3L) points in the first Brillouin zone; the 3L x 3L Hermitian 
matrices A, B and C, the last diagonal {Cjij'if = Cji 6jjf Sw) and independent of k, are 

[^(k) + I3{k) - C]^,^^.,,, = -2 Jo cos(0,r -0,7) Sy, (e^^^ 6,^i+, + e~'''^ <5,,,_i) 



-D sin( 



ji'-(t>ji) Sj'j {e'^' 6i>,i+i - e '^^ 6i',i.i) 



+6 Ji cos{(f)j'i-4>ji) Si'i (z/j, 6j'j+i + z/k ^j'j-i) , (6) 

[^(k) - B{k) - C]^.^^.,p = -2Jo6fj (e*^^(5r,/+i + e-*^^(5r,/_i) + QJi6n {ly^^fj+i + K^fj-i) , 

(7) 
Cji = 2 Jo [cos(0jy+i-0j7) + cos(0j7-0jj_i)] + D [sin(0j,i+i-0j7) + sm{(f)ji-(pj^i_i)] 

-6 Ji [cos(0j+i,/-0j7) + cos{(/)j-i^i-(t)ji)] + 18 Jihcos(f)ji , (8) 

where z/^ is the in-plane structure factor 



^'k = I {exp(iA;^) + exp ^i y-k^ + ^/Skyj + exp ^i (^-kr,-V3kyj | . 



(9) 



After a standard transformationEJ'O to creation and annihilation operators 7 ; and 7-; for 
the spin-wave excitations, 7^2 takes the form 

S N 



^2 = --:^j:c,i+sj:j:^a^) 

jl k jl 



7i(k)7,,(k) + - 



(10) 



The ground-state value is E2 = (0|7i2|0); the excitation energies e (> 0) are found from 



det 



{A -B){A + B)-eH 



. 



(11) 



The eigenvalue problem for e^;(k) is Hermitian, but the result for ej;(k) can be imaginary 
(as when the classical, in-plane solutions become unstable at intermediate fields}. 

The LSW analysis of quantum fluctuations in the C state is straightforwardE^; of course 
the antisymmetric DM term does not appear, and so quantum selection occurs just as in the 
TAFMH. For the IC state (the 111 state of Ref. ^5]), the total energy Ed + £"2 was found as 
a function of the period L and minimized with respect to L to give the optimal wavenumber 
q = 27t/L. In the relevant field region, the optimal L ranged from ^ 70 a-b plane spacings to 
~ 150, large enough to justify our treatment of the IC state as a high-order commensurate 
state and small enough that the diagonalization was practicable. 

The IC wavenumber (Figure 1) may flatten out with field, but the assumption of in-plane 
spins fails for h ^ 0.38 . More seriously, LSW theory gives a transition to the C state at 
lower field, h ^ 0.32 . The result for rj = 8.6 x 10^^ would not be qualitatively different. 

LSW theory fails because it finds the C-state energy more accurately than the IC-state 
energy. The C-state energy is found by choosing the classical configuration (for example, 
the colinear state a.t h = 1/3) which minimizes the total energy Ed + E2. Such a choice 
is not possible for the IC state and so LSW theory does not take into account sufficiently 
the breaking of the continuous degeneracy in determining the spin structure; moreover, the 
classical IC state is a poor approximation to the quantum IC state, as we now discuss. Much 
of the classical IC structureB^ can be described as a sequence of degenerate commensurate 



6 



states (a "spatially varying commensurate state"); domain walls (where the spins assume 
orientations classically forbidden in the TAFM) are well defined only very close to the IC- 
C transition. But quantum fiuctuations destroy the continuous degeneracy and so the IC 
state at even small fields will have the conventional structure, namely distinct regions where 
the order parameter is nearly constant at a commensurate value, separated by walls where 
it varies rapidly. This structure is very different from the classical structure, and thus is 
difficult to obtain by the perturbative approach of the 1/S expansion. 

Extension of the analysis to higher order (including spin rearrangement due to quantum 
fiuctuations) would also not find the two energies to the same accuracy; and it is computa- 
tionally out of reach, entirely without appeal, and even not possible if r/ = (the classical 
ground state cannot be found for h ^ 0.38). A very different approach is required. 

III. PHENOMENOLOGICAL TREATMENT OF QUANTUM FLUCTUATIONS 

The most important effect of quantum fluctuations in this problem is the breaking of 
the TAFM classical degeneracy. The following describes an innovative phenomenological 
method to treat quantum fluctuations in the TAFM part of the Hamiltonian. We represent 
the quantum part of the TAFM energy (terms proportional to S^, S^, etc.) by the term 

-E'fiuct — ~'^Z^Z^ i^in ■ Sj„) (12a) 

(ij) n 

with J2 > 0; for in-plane spins, 

jY 3 L 
-E^fluct = -J2 S^— Y^ Yl cos^{4>j+i,i-4>3i) ■ (12b) 

^ i=i 1=1 

The motivation for this form of -Efluct is partly that it favors the colinear state. A biquadratic 
coupling term is important in the classical theoryEa of the plateauE^I in the magnetization 
of the S = 7/2 system CgEu , and Ref. |^ pointed out that such a term in the Hamiltonian 
(with S > 1/2) may have the same consequences as quantum fluctuations. Our treatment 
is phenomenological because a biquadratic coupling term is not allowed in the Hamiltonian 
for CsCuCls (because 5* = 1/2), and therefore the term appears here in a fundamentally 
different way than previously. 

To test whether the phenomenological term describes qualitatively the effects of quantum 
fluctuations, we applied the formalism to the TAFM, minimizing the total energy 

3 
E = AT ^ [2 Ji S, ■ S,+i - (7 /iB H ■ S,/3 - J2 (S, • S,+i)' /S^ ] (13) 

i=i 

(this is just Eel + -Efiuct) for the flve TAFM states shown in Figure 2; the direction of H is 

irrelevant and all states have the same classical energy. 

(a) The umbrella state: The angle between the spins and the fleld is found from 

18 Ji S cos — 9 J2 S* cos (3 cos^ (p — 1) = g /xb H . (14) 



(b) The coplanar state 0i = vr, 03 = —02 with sin 02 7^ 0: The angle 02 is found from 

6 JiS (2 cos 02 — 1) — 6 J2 5* cos 02 (1 + 2 cos 202) = g ji^H . (15) 

(c) The cohnear state 0i = vr, 02 = 03 = 0. 

(d) The coplanar state 03 = 02 with sin 02 7^ 0: The angles 0i and 02 are determined by 

12 Ji S'sin(0i — 02) — 6 J2 5'sin(20i — 202) = gfisH sin (j)i , (16a) 

sin 01 + 2 sin 02 = . (16b) 

(e) The coplanar state 0i = 0, 03 = —02 with sin 02 7^ 0: The angle 02 is found from 

6 Ji5(2cos02 + 1) -6 J2 5 cos 02 (1 + 2 cos 202) = gfiBH . (17) 

State (b) is optimal at low fields, < H < Hi = 6{Ji — 3J2)S/gfiB , the region over 
which it exists. The colinear state (c) is optimal at intermediate fields, Hi < H < H2 = 
6{Ji + J2)S/gfiB', the magnetization in this region is one-third the saturation value. State 
(d) is optimal at higher fields, H2 < H < Hq, = lS{Ji — J2)S/gfiB, the region over which 
it exists. The aligned state 0j = is optimal for H > Hq . Except for the definitions of 
the fields Hi , H2 and i^s , these are just the 1/5* resulta^l for quantum fluctuations: the 
colinear state is stable in a region about h = 1/3 , states (a) and (e) are optimal at no field 
H > and < Ha, , etc. However, quantum fluctuations do not renormalizeo Ha, , and the 
expressions for the fields Hi and H2 differ qualitatively. 

Having established that the term -Efluct provides a reasonable description of quantum 
fluctuations in the TAFM, we next compare it quantitatively with the LSW energy E2 for 
the five commensurate states of CsCuCl3 derived from the TAFM states of Figure 2. The 
energies E2 were found by standard methodsEl with parameter values Jq = 28K, rj = 0, 
Ji = 4.9 K and S = 1/2. The energies -Efluct , on the other hand, are the same as for 
the TAFM and are independent of Jq . To make a fair comparison {E2 is only the leading 
quantum correction), we found -Efluct to only first order in J2 ; this means, for example, that 
COS0 = h = gfi^H/lSJiS for the umbrella state. Figure 3 plots the two energies for the 
first four states relative to state (e) (which exists over the entire field range < H < H^ 
and so is a convenient reference state). Good agreement is obtained for J2 = 0.2 K and so 
we expect that the phenomenological term -Efluct captures well enough the breaking of the 
classical degeneracy by quantum fluctuations. Because E2 is only the first-order correction, 
a moderate adjustment of J2 is acceptable; an effective J2 smaller than 0.2 K is suggested 
by the next corrections at fields just below Ha, . A separate comparison should really be 
made for rj = 8.6 x 10~^ but this seems unwarranted. 

The phenomenological approach was then applied to the IC phase of CsCuCl3 in trans- 
verse field. The total energy is E^ + £'fluct ■ The Euler-Lagrange equations are 

- 2 Jo (S,, ,_i + S,,+i) + 27] Jo (41i + 4I1) z - D (S,, ,+1 - S,,_i) X z + 6 Ji (S,-_i,z + S,+i,0 

-gfiBHk-2 Xji Sji - 6J2 S'^ [(Sj_i,i ■ Sji) Sj-i^i + (Sj+i,; ■ Sji) Sj+i^i] = ; (18a) 
for in-plane spins, these simplify to 

- 2 Jo 5 [sin(0j7-0j,,_i) + sin(0j7-0j,i+i)] - DS [cos(0j7-0j,,+i) - cos{(j)ji-(j)j^i^i)] 



+ 6JiS [sin(0j7-0j_i,i) + sm{(j)ji-(j)j+i^i)] - g ^.^H sin (pji 

-3J2S [sm(2(/.,,-20,_i,/) + sin(20,,-20,+i,/)] = . (18b) 

As described in Section II. B, scores of solutions of these equations were generated and 
followed in field, in steps Ah = 0.01 [with h = gfiB/{l8JiS)]. In the following, we use the 
value J2 = 0.13K which gives a plateau in q; there is no plateau for significantly smaller 
values, and q is nonmonotonic for significantly larger values ( J2 = 0.17 K for example)^. The 
phenomenological term does not change the zero-field spin configuration, a simple helical 
structure with wavenumber go = arctan(Z}/2Jo) ; the result for weak fields is almost the 
same as for J2 = 0. But quantum fluctuations, represented phenomenologically, make major 
differences for H ^ Hq,/3 . 

For isotropic exchange [r] = 0), the ground state is incommensurate up to Hs . The spins 
lie in the planes until h = 0.38 ; above this field, two out-of-plane solutions compete, their 
energies crossing several times; quantum fluctuations destroy the glassy phase. 

For 1] = 8.6 X 10^^, the low-field, in-plane solution is optimal for h < 0.42 ; as for the 
classical analysis (Section II. B), the weak, easy-plane anisotropy opens a large window for 
the commensurate state, which is optimal for h > 0.44 up to H = Hs ; an out-of-plane IC 
solution is optimal over such a small field range {Ah < 0.013 about h = 0.43) that we ignore 
it. Figure 4 compares the experimental results (available only to H = 13T or /i = 0.43) 
for the IC wavenumber with the theoretical values. The agreement with experiment is 
moderately good. The plateau occurs at about the observed value of q/qo and over about 
the correct field range; of course the value J2 = 0.13 K was selected to give a plateau, but 
the position of the plateau is not adjustable. The major difference is that the theoretical 
value drops prematurely. The magnetization is moderately rounded for h ^ 1/3 , but does 
not have the plateau seen in experimentJl^. 
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FIGURES 

FIG. 1. Dependence of the reduced wavenumber q/q^ on the reduced field h. The soUd circles 
are the experimental results of Ref. 20. Both the classical result (solid line, from Ref. 24) and the 
linear spin-wave result (diamonds) were found for in-plane spins and isotropic exchange (r/ = 0). 

FIG. 2. The five states, the umbrella state (a) and the four coplanar states (b) to (e), used in 
the examination of the phenomenological energy E'fluct • 

FIG. 3. Comparison of the linear spin-wave result E2 (squares) and the phenomenological 
energy -Efiuct (lines, with J2 = 0.13 K, to first order in J2) for the five commensurate CsCuCls 
states corresponding to the states of Figure 2. The figure gives the energies Eg, — E^, etc, of states 
(a) to (d) relative to that of state (e); the energy of state (b) is plotted only for h < 1/3, that of 
(c) only for h = 1/3, and that of (d) only for h > 1/3. 

FIG. 4. Dependence of the reduced wavenumber q/qo on the reduced field h. The circles are 
the experimental results of Ref. 20; the structure is unknown beyond h ~ 0.44 . The line gives the 
theoretical results based on a phenomenological treatment of quantum fluctuations (J2 = 0.13 K), 
for weak easy-plane anisotropy (r/ = 8.6 x 10~^); the commensurate state is stable for h > 0.44. 
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